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Neutrinos that decay leave their imprint on the cosmic microwave background. We calculate 
the CMB anisotropy for the full range of decaying neutrino parameter space, and investigate the 
ability of future experiments like MAP and Planck to probe decaying neutrino physics. We adopt 
two approaches: distinguishing decaying neutrino models from fiducial ACDM, and measuring neu- 
trino parameters. With temperature data alone, MAP can distinguish stable neutrino models from 
ACDM if the neutrino mass > 2 eV. Adding polarization data, > 0.5 eV is distinguishable. 
Planck can distinguish > 0.5 eV with temperature alone, and rrih > 0.25 eV with polariza- 
tion. MAP without polarization can distinguish out-of-equilibrium, early-decaying models as long 
as (mh/MeV) 2 td/sec > 230, and with polarization if (mh/MeV) 2 td/sec •> 150. For Planck with- 
out polarization, models with (m^/MeV) 2 td/sec > 9 are distinguishable, and with polarization 
if (mh/MeV) 2 td/sec > 6. Models in which neutrinos decay in equilibrium are indistinguishable 
from ACDM. Late-decaying models (10 13 sec <t d <4x 10 17 sec) are distinguishable from ACDM if 
• m h }t 5 eV for MAP and mh ^ 2 eV for Planck. Adding decaying neutrino parameters to the set of 

cosmic parameters, we calculate the statistical uncertainty in the full set of cosmic parameters. The 
ability to measure neutrino parameters depends sensitively on the decaying neutrino model. Adding 
neutrino parameters degrades the sensitivity to non-neutrino parameters; the relative amount of sen- 
sitivity degradation depends on the decaying neutrino model, but tends to decrease with increasing 
experimental sensitivity. 
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0^ ■ The anisotropy in the cosmic microwave background (CMB) can be a powerful probe of the early universe. Currently 
available data has already been used to place interesting constraints on cosmic parameters Pl-PI H> 1> H H' BL1[ II" 
[To| , [ll], [l2| , and with the advent of exquisitely sensitive satellite-based experiments like MAP [ b|[ and Planck |^4[ , it 
is possible to envision using the CMB to go beyond standard parameter estimation. Many such examples have been 
considered: detecting finite-temperature QED effects , constraining variations in the fine-structure constant jl6| , 
placing limits on lepton asymmetry |L7||, and constraining Brans-Dicke theories |18 . Another possibility is to use the 
\ CMB to probe decaying neutrinos. 

Decaying neutrinos have been considered in several cosmological contexts such as big-bang nucleosynthesis JTq| , 
|0|, ||, ||, large-scale structure formation Q ||[ fH |7| || and the CMB. The CMB anisotropy for models 



I. INTRODUCTION 



in which the neutrino decays before recombination, td t rec ~ 10 13 sec, have been calculated [^7j |2j| |30|, [H] 
Current CMB data have been used to study to late-decaying models, where td > 10 13 sec, with the result that masses 



rrih > 100 eV are mostly excluded |32j, p3fl . Late-decaying neutrinos have also been studied in the context of future 
CMB experiments M|. However, as pointed out in reference j3^], previous calculations all treated the decay radiation 
perturbations as equivalent to those of the massless neutrinos. This approximation is only valid for early-decaying 
scenarios. A systematic study of the CMB anisotropy in decaying neutrino models is needed. 

This work explores the use of future CMB observations, like the MAP and Planck experiments, to probe decaying 
neutrino physics over a large range of neutrino parameter space. It is organized as follows: First, we briefly discuss 



models of neutrino decay in Sec. II. We describe the extra steps required to calculate the CMB spectra in Sec. III. In 



Sec. IV, we briefly review cosmic parameter estimation and ruling out models in the linear regime; the formalism used 
to rule out models is developed in the Appendix. We discuss how the physics of CMB anisotropy varies as a function 
of neutrino mass and lifetime in Sec [v[ This is used to break the neutrino parameter space into regions where the 
physics of the CMB anisotropy is similar. We then present results: distinguishing decaying neutrino models from 
standard models, and measuring cosmic parameters in Sec. |vi| . 



' r-lopez@uchicago.edu ; 



this paper represents partial fulfillment of the requirements for the Ph.D. in Physics at the University of Chicago 



2 



II. MASSIVE NEUTRINOS 

The evidence for neutrino mass from atmospheric, solar and direct-beam neutrino oscillation experiments is 
compelling [^5), and massive neutrinos tend to decay unless protected by some symmetry. It is therefore interest- 
ing to consider the cosmological signature of decaying neutrinos. 

In this work we will consider neutrinos that decay non-radiatively into light decay products. By non-radiative we 
mean that the decay products are electromagnetically non-interacting. Radiative decay channels could also exist, e.g., 
Vh — * or vh ~^ e + e~i>i. However, these models are generally excluded by observations unless the lifetimes are 
extremely long |35| [}(|; the region of parameter space that can be probed by the CMB is certainly excluded. There 
are several models with non-radiative decay products that are motivated by particle physics. For example, familon 
models |37], ||^, |3^] predict the following decay process: i>h — » vicj) where is a familon, a massless Nambu-Goldstone 
boson associated with spontaneous breaking of a continuous, global family symmetry. In these models, the decaying 
neutrino mass and mean lifetime are related at tree level by 

16tt / F x ' 



U = — — , (1) 
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where F is the energy scale at which the family symmetry is broken, and it is assumed that the neutrino vi is 
much lighter. This interaction induces a corresponding charged-lepton decay, and experimental constraints on their 
branching ratios can be used to set lower bounds on F . Familons corresponding to a fi-r family symmetry are the 
least well constrained: the branching ratio B(t — > yap) < 3 x 10~ 3 |4(| which implies that F>4x 10 6 GeV for the 
second-third family symmetry. This leads to the following constraint, assuming that = v r and /// = v^: 



*L) (^) 3 >3 .0xl0 17 . 
sec / V ev / 



(2) 



Much of the decaying neutrino parameter space that can be probed by the CMB satisfies this constraint. 

In models where neutrinos acquire mass through the see-saw mechanism, the three-body decay v>h — * vivi&i can 
occur |56), mediated by the exchange of a Z-boson. However, the lifetime for this decay, 

td 1 O 30 ( eV 



— -10" — , (3) 
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is so large that the neutrino is effectively stable over the interesting region of parameter space. 

Motivated by the discussion above, we consider the following decay channel throughout this work: Uh ~* vi<fi. 
However, alternate decay processes, like the aforementioned Vh — > viViV\ do not alter our results much; the small 
differences are discussed where they exist. Therefore, it is appropriate to specify decaying neutrino models by 
and td alone. The results are then model-independent for most of the interesting parameter space. 

III. CALCULATING THE ANISOTROPY 

The anisotropy in the effective temperature of the CMB radiation, ST, is typically described in terms of spherical 
harmonics, 

^I^E^cm). (4) 

Im 

where 9 and <f> describe the position on the sky, and To = 2.728 K is the mean background temperature of the CMB. A 
given theory, specified by some set of cosmic parameters, makes predictions about the distribution of the coefficients 
aj m . For Gaussian theories like inflation, the coefficients are drawn from a normal distribution, with zero mean. In 
this case, all of the predictions of the theory are encoded in their variance. Therefore, the predictions of the theory 
can be written in terms of the C; coefficients, defined by 

C Tl ^{aJ m aJ; n ) . (5) 



In general, the temperature anisotropy does not contain all of the information in the CMB because the CMB is 
polarized. The symmetric, trace-free polarization tensor Fab can be decomposed into two kinds of scalar modes with 
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opposite parities: an electric-type mode and a magnetic-type mode The polarization field can be expanded in 
terms of electric and magnetic type spherical harmonics 5^ m V^ 6 \ 5 with parity (—1)' and (— respectively: 



Tn 



Eh 

Irn 



FmYlm(ab)(0' ' 



a lmXlm(ab)(^' ' 



(6) 



When polarization is included, the information in the CMB anisotropy can be characterized by three additional 
correlation functions 
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Because the magnetic mode has parity opposite the temperature and electric modes, the T — B and E — B correlation 
functions vanish [ [42] . In this work we assume that the primordial perturbations are purely scalar density perturbations, 
with no tensor component. Their lack of handedness implies that scalar density perturbations cannot generate the 
magnetic- type modes [43 . Therefore Cbi = for the models we will consider. This assumption is motivated by 
the fact that most inflationary models produce tensor fluctuations too small to be easily detected, even with future 
satellite-based experiments Q|. In any case, for simplicity we will ignore this possibility. 

The CMB anisotropy is related to perturbations to the photon distribution function, which is itself coupled to 
other particle species and gravitational metric perturbations through particle interactions and gravity. In this work 
we use the synchronous-gauge, where the coordinate and proper time of freely-falling observers coincide; all of the 
metric fluctuations occur in the spatial part of the metric, ds 2 = a 2 (r) [— dr 2 + (Sij + hij) dx % dxi~\ . The metric 
perturbations can be decomposed into scalar, vector and tensor components; we will be concerned solely with the 
scalar perturbations. These can be written in terms of two scalar functions h and 77 |}45|. In Fourier space, 



hij(k, t) — |fcifcj h(k, t) 



(8) 



where k is the Fourier mode and r is conformal time defined in terms of regular time t and the scale factor a by the 
relation dr = dt/a. To calculate the CMB anisotropy we need to know the metric perturbations h and rj, as well as 
the distribution functions for all components: decaying neutrinos, decay radiation, photons, massless neutrinos and 
cold dark matter (CDM). The differences between a standard scenario with no decaying neutrinos, and the decaying 
neutrino scenarios we consider can be summarized as follows. In a decaying neutrino model: 

• The energy densities of some of the components evolve differently from the standard case. This affects the 
dynamics of the expansion of the universe thr ough t he Fricdm ann equation, i.e., the Hubble parameter a/ a is 
modified. This modification is covered in Sees. Ill A and III B . 



• The Boltzmann equations that govern the evolution of the decayi ng ne utrino and decay radiation perturbations 
must be modified to include decay terms. This is covered in Sec. Ill C . 



A. Friedmann equation 



The evolution of the scale factor is governed by the Friedmann equation. For the flat universes considered here, 

(a) (a ) 3m|^ a ^' ^ 

where m p = 1.221 x 10 22 MeV is the Plank mass and p{a) is the total energy density. In this work, overdots are used 
to denote derivatives with respect to conformal time. The total density can be broken into components: the decaying 
neutrino ph, its decay products p r d, standard radiation p sr , i.e., photons and two massless species of neutrinos, CDM 
+ baryons p m , and vacuum energy density p\. The standard components evolve simply with scale factor: p sr oc a~ 4 , 
p m oc a -3 , p\ oc a . However, decays (and possible inverse decays) complicate the decaying neutrino and decay 
product density evolution, which complicates the Friedmann equation and makes it impossible to solve analytically, 
except in special cases. 
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B. Energy density evolution equations 



The distribution function for the i-th component, fi (x J , q J , t) depends in general on seven variables: position x 3 , 
comoving momentum q\ — ap 3 where pi is the proper momentum, and conformal time r; it evolves according to the 
Boltzmann equation, 

dfi _ dfi dx j dfi dn j &U dq dfi _ 

where n 3 is a normalized vector in the direction of the momentum, q 3 = qin 3 , and C[/j°] is a collision functional that 
describes particle interactions. The factor of a multiplying the collision functional is just convention; it is a conversion 
between conformal time and real time, where collision terms are more easily described. 

To find the equations governing the evolution of the energy densities, we consider the Boltzmann equation for the 
zeroth order distribution function, ff{q, t), denoted with a superscript-O. By zeroth order, we mean that we are 
neglecting the spatial perturbations in the distribution functions, so that the term proportional to dff jdx 1 = 0. The 
quantity dq/dr is first order in the metric perturbations |^5|, so that it too can be neglected. We also assume that 
/j° does not depend on the momentum direction (df^/dn 3 = 0), but allow to have arbitrary dependence on qi. In 
this limit the Boltzmann equation simplifies: 

-£r=aC[f?]. (11) 

For the decay process Vh — * v\ <f), the component i is either the decaying neutrino (i — > h) or one of the decay products 
(i — * I, or i — * <j>). We can find the zeroth-order energy density from the distribution function using the definition 

= d ^^ ] - < 12 ' 



where = y/qf + a 2 mf and <?i is the number of internal degrees of freedom for particle i. For the massive and 
massless neutrinos, = gi = 2, and for the scalar decay particle <ju = 1 since it is assumed to be spin-0 and its own 
antiparticle. 

We next turn to the collision terms. In general, every type of interaction that the particle experiences will contribute 
to the these terms. Fortunately, in the case of decaying neutrinos, only a few interactions are important. Because 
the decaying neutrino interacts with the rest of the universe via the weak interaction, it decouples at a very high 
temperature of order a few MeV, just like standard, massless neutrinos. So for temperatures of interest here (eV-scale 
rather than MeV-scale), the decaying neutrino-decay radiation system is decoupled from the rest of the universe. 
Therefore, the only processes that are important are decays and inverse decays. Scatterings can be neglected in the 
calculation of energy densities, since they just shuffle energy among particl es [|lT 

For the massive neutrino (i — * h) the collision functional can be written |H 



C h [f h ] = -V D +Y 
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r h D = i^AV) / d Ql [i + fl{ qi )] [i - f?(e h - qi )] , 

T h iD = r— C 1 - fnilh)] / d gi $(qi)f?(e h qi) . (13) 

tde h q h Jl/2(e h -q h ) 

In this expression T D arises from decays, — > vi(f> and Tj D arises from inverse decays, vi<f> — > v^. The integration 
limits follow from the kinematics of the interactions. 

The collision terms for the decay products are similar to those for the decaying neutrino. For the light neutrino 

(i-Z), 



Ci[f?] = -r l D + r l ID , 
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(14) 
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and, for the scalar particle (i 



CM 



VI = 2 
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The Boltzmann equation for each type of particle, Eq. O, their collision term equations, Eqns. |T^-|T^, and the 
Freidman equation, Eq. |^, determine the dynamics of the expansion of the universe. They form a closed set of 
integro-differential equations for the evolution of the scale factor, and require numerical methods for their solution. 
In particular, the collision term integrals are complicated functions of momentum. However, for certain special cases 
these equations simplify, and for other cases we can estimate the late-time densities without having to solve the 
equations at all. 



1. Out-of-equilibnum decays 

Neutrino decays become important when the age of the universe is near the neutrino mean lifetime. If T(td) <C 
where T(td) is the temperature of the universe at time td after the big-bang, the neutrino decays non-relativistically, 
so that when the neutrino starts to decay, the thermal energy of the decay products cannot overcome the rest mass 
energy of the decaying neutrinos. This suppresses inverse decays relative to decays and causes the decays to occur out 
of equilibrium. We will use the terms out-of-equilibrium decays and non-relativistic decays interchangeably. Thus, 
the neutrino decays away when t ~ td- These decays can generate a large amount of decay radiation, depending on 
the initial abundance of the decaying neutrino and how non-relativistic the neutrino is at decay. 

Since neutrinos decouple from the rest of the universe at a very high temperature determined by their weak 
interactions, all of the neutrinos, including the massive, decaying neutrino, are ultra-relativistic at decoupling (we will 
not consider MeV-scale decaying neutrinos). Thus, their abundances are large, of order the photon abundance. The 
decaying neutrino, if still present, becomes non-relativistic when T < Then its energy density scales as matter, 

as a~ 3 instead of as radiation, which scales as a~ 4 . Its energy density, and consequently the energy density of its decay 
products, becomes relatively more important the longer the decaying neutrino is still around and non-relativistic. 

For out-of-equilibrium decays, simplified evolution equations for the decaying neutrino and decay radiation densities 
can be found. The collision term for the decaying neutrino simplifies, because in this limit we can neglect ff and /°: 

Ch [f h ]^^f h . (16) 

The Boltzmann equation can then be converted into a differential equation for p^, by multiplying each term by p\Eh 
and integrating out ph- We find that 

d -f + Z a -( P UP° h ) = -^nl, (17) 
or a v td 

where P° is the pressure and n° is the number density of the decaying neutrinos, given by the definitions, 



1 a 4 2ir 2 
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/ dq^-f^q). (18) 



A couple of comments about the evolution equation for pi are in order. Note the presence of a pressure term on 
the left hand side. In the limit of completely non-relativistic decays, this term is zero, but otherwise this term can 
be a significant correction. If we neglect the pressure term, then the 3(a/a)ph term represents the fact that matter 
density varies as a -3 in the absence of decays. A similar comment applies to the product mh on the left hand side. 
For completely non-relativistic decays, all of the decaying neutrino energy is rest mass energy so that m^nl = p°, 
but otherwise the two quantities are not equal. 
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Given the decaying neutrino energy density, the decay product energy density p rd = pi + p^ can be obtained from 
the first law of thermodynamics : 

d(a 3 p rd ) u da 3 d(a 3 p h ) 

-^^-- Prd ^7~^^- (19) 

Since the decay radiation is massless, P rd = l/3p r d, and we find that 

dprd .,a ( dp h a \ a 

— h &~Prd = — I —j 1" 3-ph 1 = — —Ph > (20) 

or a \ cLt a J td 

where the second equality holds for fully non-relativistic decays. In the absence of decays, this equation implies that 
the decay radiation density scales as a~ 4 , as expected for massless particles. Finally, we can obtain a simpler equation, 
which will be useful later, for the evolution of the decay radiation density. Let r r d = Prd/Piv where p\ v is the density 
in a single species of standard, massless neutrinos. Then Eq. ^0|, and the fact that p\ v oc a -4 , implies that 

dr rd m h n° h a 

-j— = T' ( 21 ) 

ar p\ v td 

To find the energy densities of the decaying neutrino and its decay radiation for out-of-equilibrium decays, we 
numerically solve Eqns. |lj and |2(], together with the Friedmann equation, Eq. || ^6|. Results for several decaying 
neutrino models are shown in Fig's [j] and ^|. There we plot the energy densities, scaled by the critical density, for 
all of the components: standard radiation, CDM, vacuum energy density, decaying neutrino and decay radiation. 
The first figure shows a succession of masses with lifetimes fixed at 10 9 sec. These are models where the neutrino 
decays before last scattering, t rec ~ I0 13 sec. It is easy to see that the decay radiation becomes more important as 
the mass increases, in keeping with Eq. ^3|. If the neutrino is massive enough, then it can cause an early phase of 
matter domination before it decays and its decay radiation dominates. The second figure shows some models where 
the neutrino decays after last scattering. 

By determining how non-relativistic the neutrino is when it decays, it is possible to obtain an estimate of the 
energy density in decay radiation, without resorting to the full Boltzmann equations. To do this we define a relativity 
parameter a, proportional to the square of the decaying neutrino's mass divided by its thermal energy at the time 
of decay, and with the property that a ~ 1 at the border between relativistic and non-relativistic decays. For non- 
relativistic decays, a ^> I and p rd / 1 p\ v is large; for ultra-relativistic decays, a -C I and p rd / 1 p\ v ~ f . Consider a scenario 
with a ~ f . Here, the universe is never dominated by the massive neutrino. For a radiation dominated universe at 
decay, the Friedmann equation gives the relation between decay time and temperature Jj^l , t d — 0.3(7* ^"W^J; 
where ~ 3.36 is the effective number of relativistic degrees of freedom. Since the temperature at decay Td ~ to/j/3, 
the neutrino parameters enter in the combination m\t d , which implies that 



0.11 



f"*YlU\ (22) 
VMeV/ Vsec/ v ' 



Because matter density decreases as one power of the scale factor relative to radiation density, we can estimate the 
energy density in decay radiation in units of standard massless neutrinos, N rd , as follows: N rd — 0,4/ a nr , where a nr 
is the scale factor when the neutrino becomes non-relativistic and ad is the scale factor at decay. Here we assume 
that the decay instantaneously transforms the density in decaying neutrinos to the decay radiation. If the universe is 
dominated by the decaying neutrino at decay, then the Friedmann equation can be used to obtain a d . The result is 
that Q 

N rd ~ 0.52 a 2/3 , (23) 



valid for a ^> 1. The numerical coefficients in Eqns. [22| and 23, but not the overall dependence, have been fudged 
by a small amount so that the the formula for N r d agrees well with numerical results. The bottom pane of Fig. || 
shows this estimate versus numerical results for the total radiation density N v — N rd + 2, as a function of a. The 2 
represents the two species of massless neutrinos. As the figure shows, the agreement is good as a rough estimate. 



2. Equilibrium decays 



If T{td) S> rrih/3, neutrino decays become important while the neutrino is still ultra-relativistic. In this case both 
decays and inverse decays occur, and the collision terms do not simplify. It is, however, possible to obtain an estimate 
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t d = le+09, m h = le+02 t d = le+09, m h = le+03 




log a log a 

FIG. 1: The evolution of the energy densities, relative to the critical density, of the various components of the universe, in early- 
decaying scenarios. The notation is as follows: solid line = CDM+baryons, long-dashed line = standard radiation (photons + 
2 massless neutrinos), dotted line = decaying neutrino, dot-dashed line = decay radiation. For all values of the scale factor, 
53 ; Qi = 1. The background cosmological model has a cosmological constant Q,a ~ 0.7 today. The vertical line represents the 
epoch of recombination. The models shown here all have td = 10 sec; mh varies from 10 2 to 10 s eV. For ran = 10 2 eV, the 
decay is barely non-relativistic: a — 1.1. The decay radiation density never matches the density in standard radiation. For 
the higher-mass scenarios, the decays are out-of-equilibrium and the decay radiation dominates the standard radiation for all 
times after decay. Another feature to be noted is the relative importance of components at recombination; this determines the 
amount of the early-ISW effect. For = 10 4 , 10 s eV, the universe is radiation dominated at last scattering, creating a large 
early-ISW effect. 



for the energy density in decay products long after the neutrino has decayed away. This estimate relies on the fact 
that when t > td,, the decay and inverse decay processes are sufficiently fast relative to the expansion rate to establish 
chemical equilibrium between the decaying neutrino and its decay products [ p2| . Then the distribution functions are 
approximately thermal in form with pseudo-temperature T 1 not necessarily equal to the temperature of the universe. 
Therefore, we have 

f° = e (Bi- w )/r' ± i ' ( 24 ) 
for i = h, I, or cj>, with pseudo-chemical potentials related by 

Mft = W + ■ (25) 
It is easy to show that the Boltzmann equations imply that the following relations hold generally : 

d 



(It 
d_ 



(a 3 (n h + nij) = , (26) 
(a 3 (n h + rty)) = , (27) 



and that for T ^> mh, the following holds: 



d 



- r (a i (p h + pi + pt))=0. (28) 



dr 
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t ri = le+15, m h = 5e+00 



t H = le+15, m h = 5e+01 
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FIG. 2: Same as Fig. [l], but for late-decaying neutrinos. In neither mj, = 5 eV scenario does either the massive neutrino or 
its decay products ever dominate the energy density. In both scenarios with mh = 50 eV, the universe is dominated by the 
massive neutrino at recombination, and by the decay radiation at decay. 




FIG. 3: Late-time asymptotic behavior of extra radiation energy density expressed in units of species of massless neutrinos, as 
a function of relativity parameter a. Plotted are the non-equilibrium limit SN V = 0.52a 3//2 , valid for a > 1, as well as data for 
2-body decays (i>h — » vi <f)) and 3-body decays (yu — > vivivi). In both the 2 and 3-body decay scenarios we have assumed an 
initial thermal abundance of heavy and light neutrinos at the standard neutrino temperature T v = (4/ll) 1//3 T 7 . 
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This equation implies that the total comoving energy density in the decaying neutrino-decay radiation system is 
unchanged by the decay /inverse decay processes, for T(i^) rrih- 

We numerically solve Eqns. |25|, 26, and ^ for fih, fii, pb$, and T'. For initial conditions we assume a thermal 



initial abundance of heavy and light neutrinos, and no initial scalar particles. We find that 

Hh = 0.092 T v , m = 0.581 T v , ^ = -0.489 T v , T' = 0.884 T v , (29) 

where T v = (4/1 f) 1 / 3 T is the standard neutrino temperature. This solution if valid for t>td and < T(i). When 
the universe cools enough so that T < JW/,/3, the inverse decays become suppressed and decays predominate. The 
rest mass of the heavy neutrinos is starting to become important, increasing their total energy relative to the massless 
case. As they decay away, the energy density in heavy neutrinos is then transferred to the decay radiation, raising its 
temperature. We can calculate the amount of heating by using the fact that the entropy of the heavy neutrino-decay 
radiation system is conserved. We find that the neutrino decays raise the decay radiation temperature by 14.7%. Our 
final result, the energy density in decay radiation, can be expressed in units of standard, massless neutrino energy 
density: N r d — 2.17. In a scenario with no decaying neutrinos, this number would be 2.0, so 8N V = 0.17, where 
8N„ is the change in radiation density. A similar procedure could also be repeated for the case of three body decays: 
Vh — > 3 2/;. In this case, 8N V = 0.52 for a <C 1, and the large-a behavior of the radiation density is identical to the 
two-body case. 

Results for the decay radiation energy density, for both equilibrium and out-of-equilibrium decays, are shown in 
Fig. |^. To summarize, for T(td) <SC rrih/3, decays occur in equilibrium, with decays and inverse decays important for 
t^td- For td < t < t(T = rn/j/3) the energy density in radiation is repartitioned, but the total value is the same as 
if the neutrino did not decay. For t > t(T = m^/3), the heavy neutrino decays away and increases the total radiation 
density by 2.3%. 



C. Perturbation Boltzmann equations 



Following reference p5| , we would like to derive a hierarchy of Boltzmann equations describing the evolution of 
perturbations to the decaying neutrinos and the decay radiation. The i-th distribution function can be written as the 
product of an unperturbed, thermal function times a perturbation, as follows: 



fi (ar* , q, n> , r) = /»(«,, r) [l ± ^ (x* , q, n* , r)] . 



(30) 



The Boltzmann equation, Eq. |lC|, then yields an equation for the evolution of the perturbation. Upon taking the 
Fourier transform, 
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(31) 



Since the decay radiation is effectively massless, e r d = q r d, and we can integrate the momentum dependence out of 
the Boltzmann equation. We define a momentum-independent perturbation F r d as in reference Ej| , scaling it by the 
decay related factor r r d for convenience: 



,r , = Jdqq 3 f? d (q,T)y rd (k,q,n,T) 
r r d\h,, it, I — . Q I r d ■ 

J dqq 3 f» d (q,T) 



(32) 



Unfortunately, the complicated form for the collision terms in the Boltzmann equation makes it difficult to derive 
simple equations in the general case. For the rest of this section, we will specialize to the case of out-of-equilibrium 
decays, where these terms simplify. Then the Boltzmann equation governing F r d can be shown to be p2| 



F rd + ikfj,F rd + 4 



h + 6f) 



Trd 



= r rd N , 



(33) 



where 



/ dq h qlfhilh, T)^ h (k, q h ,r) 



N {k,T) = 



8 / q h 



Jdq h q 2 h f^(q h ,r) 



(34) 
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fi = k -h and P n {n) are the Legendre polynomials of order n. In Eqns. |3^ and |34|, only terms up to 0(qf i /a 2 mf l ) have 
been kept. Similar equations for the evolution of perturbations in the decay radiation can be found in references [p6| , 
©and g. 

The dependence of F r d on \x can be eliminated by expressing it as a series of Legendre polynomials, F rc i = F r dd Pi, 
leading to the following Boltzmann hierarchy for the decaying neutrino perturbations, valid for out-of-equilibrium 
decays: 

Kd = — \ (h + 28 r d) - — (3rd — Sh) , 

3 V / r rd 

@rd = k 2 ( -5 r d — &rd) —Ord j 

V 4 J r rd 

2 / „ ; _.\ r 
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26 rd + h + 6rj) - —<J r d 
> r rd 



rd.l 



[lFTd,i-x-{l + l)Frd t i + i]-—F rdtU l>3 (35) 



21+1 T rd 

where S r d — F r dfl/r r d, Q r d = 3kF r d,i/4r r d and a r d = F r d.2/2r r d- This set of equations is identical to the Boltzmann 
hierarchy for standard massless neutrinos, with the addition of the terms proportional to r r d/r r d oc l/t d . The extra 
decay term can have a large effect on the decay radiation perturbations when t ~ t<j> but f° r late times the perturbations 
approach the values they would have attained in its absence. To calculate the decay radiation perturbations, we added 
a separate Boltzmann hierarchy, described by Eq. In our numerical scheme, this hierarchy must be terminated at 
some value of multipole moment l en d- We terminate the hierarchy by adding the extra equation for F r d,i snd +x, 

2l en d ~t" 1 

F r d,i end +i = — -j— — F rdi i eni - F rd ,i and -i , (36) 

the method used for the massless neutrino hierarchy in CMBFAST |47[ ]. 

Because the decaying neutrinos are massive, th ^ qhi and it is impossible to integrate the momentum dependence 
from their perturbations. After expanding the decaying neutrino perturbation in terms of Legendre polynomials, 
= Yli ^h,iPi, the Boltzmann equation becomes 

^h,o = + ^h— ^ h ,o , 

eh 6 dmqh t d th 

t Qhk , T . arrih 

• qhk , . / 1 ; 2 \ din ami T 

*fc,a = (2* h ,i - 3* /li3 - \—h + -j) — - -^* h ,2 , 
oe^ \15 5 / dln% f d e/, 

* W = e + ^ +1) [Wm-x - (* + D*M +1 ] - , I > 3 ■ (37) 

This set of equations differs from the evolution equations for massive, non-decaying neutrinos only through the 
presence of the term proportional to l/td- The decay term is easily interpreted. For non-relativistic neutrinos, the mh 
in the numerator cancels the eh in the denominator; the result is just the differential equation for exponential decay, 
in conformal time. If the neutrinos are not completely non-relativistic, then their velocities become important, and 
there is a time dilation factor associated with transforming between the neutrino rest frame and the thermal frame. 
In this case, mh/eh becomes the special relativistic gamma factor for this transformation. 

It should be noted that, for out-of-equilibrium decays, the perturbation evolution equations are independent of the 
details of the decay radiation, except for the fact that is must be light and weakly-interacting. The energy density 
equations, Eq. ^ and [l?], are also independent of the details. Therefore, the CMB anisotropy becomes a function 
of nih and td, independent of the decay channel. In fact, the calculations can be generalized to encompass generic 
decaying particles. The main difference in the generic scenario will be due to the initial abundance of the decaying 
particle which will depend on its interactions. However, a generic decaying particle will produce a CMB spectrum 
very similar to a decaying neutrino with the same lifetime, provided that the densities of decay radiation are the same. 
Finally, note that this simplification is valid for out-of-equilibrium decays only. 



IV. ANALYZING THE DATA 



This section briefly reviews estimating cosmic parameter uncertainties ("error forecasting"), and using data to rule 
out decaying neutrino parameter space. For further discussion of error forecasting in parameter estimation, see e.g., 
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Refs. []42|, |4g, [49], |5C], plL |5^]. The formalism used to determine which decaying neutrino models can be ruled out is 
discussed in the Appendix. 

A. Measuring m;, and td 

A given theory, specified by a set of cosmological parameters {A^} (i — 1 . . . N, with N the number of cosmic 
parameters considered) makes predictions about the multipole amplitudes, the C;'s. The results of a CMB experiment 
are estimates of the CVs, with some experimental uncertainties. Of course, we cannot know in advance the values of 
C;'s that a given experiment will measure; however, by knowing what we expect for the uncertainties, we can estimate 
how large the uncertainties in the parameters should be. 

For an experiment with data out to some maximum I — l m ax, we can define a goodness of fit statistic that is a 
function of {Aj}: 

X 2 ({AJ) = £ £ [cT rV {{^}) - C d xT] Vxyi [C?r v i{^}) ~ C d Y f a ] , (38) 

1=2 X,Y=T,E,C 

where C^"^ is the theoretical spectrum for cosmic parameters {A^}, C^ a ; ta is the measured spectrum and Vxyi is 
the covariance matrix between estimators of the different spectra. For a cosmic variance limited experiment with data 
to some maximum I = l max , the diagonal components of Vxyi are given by PS 



2 

Vtti = 27+1 Tl ' 
2 



Veei - n , ; C 2 
2 



El ! 



2Z + 1 

2 



Vcci = (Chi + CtiCei) , (39) 

and the non-zero off-diagonal components are given by 

Vtei = 2T+T C ^ 7 ' 
Vtci = 2i T ^ CtiCci , 

Veci = -^1 C eiCci , (40) 

for / lmax- 

The measured cosmic parameters, {A^}, are determined by minimizing x 2 ({Ai}): 

dx 2 



V Hi / I- 



DA, 



({AJ})=0, (41) 



for j = 1...N. If we assume that the measured cosmic parameters are close to their actual values, denoted {A^}, 
then \ 2 can be expanded about its minimum as follows: 

X 2 ({AJ) ^ X 2 ({~Xi}) + £ - ~\) Oij (A, - A,) , (42) 

where otjk is the Fisher matrix, 



I XY 



p. n theory p ir ttheory 



3 



The Fisher matrix determines how rapidly x 2 increases as the cosmic parameters are varied away from their true values. 
Under certain reasonable assumptions (53) , the uncertainties on the parameters are determined by this matrix. If we 
allow all cosmic parameters to vary simultaneously, then 



6Xi = y/ (44) 
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The formalism above assumes data for both temperature and polarization. If only temperature data is obtained, 
then the covariance matrix Vxyi becomes a number: 

2C 2 

Vxyi = 5 XT Syt , (45) 

where S is the discrete delta function. 

To calculate the uncertainties in the parameters, we will assume some decaying neutrino scenario. The set of cosmic 
parameters will include neutrino parameters, like and t^. The uncertainties will then depend on the model we 
assume and the parameters we allow to vary. 



B. Ruling out models 



It could also be the case that no theoretical model can specify the data. For instance, in a decaying neutrino scenario, 
the data could be analyzed without considering neutrino parameters. In general, two things will then happen. 1) 
The best-fit parameters will be systematically offset from the true parameters. 2) No theoretical model will fit the 
data well, i.e., the best-fit x 2 w iU be higher than expected. In special cases, one or the other thing will happen. 
For instance, if the effects on the CMB of the decaying neutrinos and their decay radiation is exactly mimicked by 
some perturbation to the set of cosmic parameters, then a ACDM model with offset parameters will fit the data well. 
If, on the other hand, the effects of the decaying neutrinos and the decay radiation are orthogonal to the effects of 
parameter offsets, then the offsets will be small, but no model will fit the data well. If no ACDM model can reproduce 
a decaying neutrino model, in the sense that the best-fit % 2 is large, then the decaying neutrino model is said to be 
distinguishable from ACDM. 



If the offsets are small, then the problem can be analyzed analytically. This is done in Appendix VII . The procedure 
we use to determine the distinguishability of a model is to calculate the Ci spectrum and the Fisher matrix for the 
cosmic parameters being considered, for the baseline ACDM model. Then, for a given decaying neutrino model we 

• Find the parameter offsets using Eq. ^2|. 

• Determine the probability distribution for the goodness of fit x 2 ■ Being approximately Gaussian, this distribution 
is fully characterized by expected the best-fit (Xmin)j given by Eq. |38| and the variance a x , given by Eq. 

• Determine the level of distinguishability by convolving the probability distribution for Xmin with the allowed 
level for each Xmiw as P er Eq- TO 



V. REGIONS OF PARAMETER SPACE 



In decaying neutrino scenarios, the physics of the neutrino decays, and therefore the CMB anisotropy changes as 
the neutrino parameters are varied. It is therefore useful to break the parameter space into regions and consider each 
region separately. To do so we first note that several physical scales naturally divide the parameter space: 

• id = tjj- This represents the division between stable and unstable neutrinos. 

• floh 2 — 0.25: For our decaying neutrino models, we let CIcdm vary to enforce a flat universe: CIq = 1. 
If the density in neutrinos or decay radiation today is large enough, then fio > 1 even with no CDM. For 
reasonable values of ho, regions with Oo > 1 tend to produce universes young enough to violate independent 
age constraints: floh 2 < 0.25. For stable neutrinos, this translates into the well-known bound on the sum of the 
masses, mi < 24 eV, where the index i runs over all neutrino species. 

• td — tree'- The decay radiation for neutrinos that decay before last scattering sources CMB anisotropy through 
early-ISW effect, while the decay radiation for neutrinos that decay after last scattering creates a late-ISW 
effect. 

• mh = 3T(t<j): This scale divides equilibrium (m/j < 3T(i<j)) from out-of-equilibrium (m/j > 3T(td)) decaying 
neutrinos. Neutrinos that decay in equilibrium produce small changes in the radiation density, while those that 
decay out-of-equilibrium produce larger effects. 

• uifi — 3T(t rec ): This scale determines whether the decaying neutrinos are relativistic (mh ^ 3T(i rec ) or non- 
relativistic (m^ > 3T(t rec ) at last scattering. 
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FIG. 4: Decaying neutrino parameter space, divided into regions according to the physics of the CMB anisotropy. 



Based on these scales, we have broken the decaying neutrino parameter space into regions, labeled alphabetically, 
as shown in Fig. |J: 

• A: 3T(t rec ) <m h < 3T(t d ) 

In this region the neutrino decays in equilibrium, before last scattering. The energy density in radiation is 
increased relative to the standard case by SN U = 0.17. The only difference between the CMB anisotropy of 
these models and the baseline model is due to this extra radiation. If the universe is not completely matter- 
dominated at last scattering, then the gravitational potentials are decaying at last scattering, when the primary 
anisotropy is being formed. Decaying potentials at last scattering generate anisotropy through the early-ISW 
effect. The small amount of extra radiation in these models induces a small amount of extra anisotropy. The 
angular scale of the effect is determined by the sound horizon at last scattering, placing the feature near the 
first acoustic peak, which, for the flat universes that we consider, I ~ 200. The degeneracy in rrih and t<j means 
that these models can be considered as a group. 

Because the CMB anisotropy in this region depends only on the radiation density at last scattering, the details 
of the decay channel are unimportant, except to the extent that they determine this density. For example, it is 
easy to extend the analysis to include the three-body decay scenario i>h — > vi v\ vi , because we know that in this 
scenario, 5N V = 0.52. 

The claim that we can calculate the CMB spectrum for models in region-A by simply adding 0.17 species of 
massless neutrinos bears examination. One possible concern follows from the fact that if massive neutrinos 
are present near last scattering, then they will affect the CMB anisotropy. However, in this region there 
are no massive neutrinos left at last scattering; they have decayed away by then. A more serious concern 
involves spatial perturbations to the decay radiation. Treating the decaying radiation by simply increasing the 
effective number of massless neutrino species effectively assumes that the decay radiation perturbations are 
equal to massless neutrino perturbations. But for times much later than those when decays are important, the 
decay radiation perturbations approach those for massless neutrinos. This is because the collision term in the 



Boltzmann equation that describes the perturbation evolution, described in Sec. QIC, is only important when 
decays are important, and the evolution equations without the collision term are identical to those for standard 
massless neutrinos. In this region, the decaying neutrinos decay away when they become non-relativistic, when 
T{t) < mh/2>. If this time is much earlier than recombination, i.e., if mh <C 3T(t rec ), then the decay radiation 
perturbations can be approximated as standard massless neutrinos, and the arguments in the last paragraph 
hold. From Fig. ^, this condition holds in region- A, for points a good deal to the right of the defining line 
n%h = 3T(t rec ). We will assume that this is true for the rest of this work. 

B: m h > 3T(t d ), t d < t rec 

Here, neutrinos decay out-of-equilibrium, before last scattering. Thus, as for region-A, the decay radiation 
sources the early-ISW effect which results in extra anisotropy near the first acoustic peak. But the effects are 
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larger in this region since out-of-equilibrium decays can generate large amounts of decay radiation, as shown in 
Eq. ^3]. The amount of extra radiation, and hence the CMB spectrum, depends on one parameter only, either a 
or SN V . This is in contrast this to the constant effect in region-A Some models from region B, parameterized 
by SN U , are shown in Fig. ||. Another effect is visible in addition to the early-ISW acoustic peak enhancement: 
a shift of all features to smaller angular scales. This is due to the fact that, as the amount of radiation at 
last scattering increases, the sound horizon at last scattering decreases. For the standard ACDM model, the 
universe is mostly matter-dominated at last scattering, with r = 2 v / eLff ~ 1 . In the limit of a completely radiation- 
dominated universe at last scattering, this relation is modified to become r = aH^ 1 . This is the reason for the 
shift to smaller angular scales, since a at last scattering is the same in both scenarios. 
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FIG. 5: CMB anisotropies for early-decaying models, corresponding to region- B of parameter space. The quadrupole-normalized 
anisotropy is plotted as a function of I. The solid line depicts ACDM. The dashed line represents a model with SN r d ~ 1.0, the 
dashed-dotted has SN r d = 10.0 and the dotted has SN r d = 100.0. Each value of SN r d corresponds to some a through Eq. ^j, 
and thus to a one-parameter family of decaying neutrino models. Note that as SN r d increase, the universe becomes less and less 
matter-dominated at recombination and the early-ISW peak becomes more prominent. For 5N r d = 100, the universe is very 
radiation-dominated at last scattering. This changes the age of the universe at recombination and shifts features to smaller 
values of I. 



In the future, we will parameterize models in region- B in terms of the decay radiation density 5N V . We should 
question the validity of this parameterization. We would expect that the complicating effect from massive 
neutrinos being present at last scattering is absent for td <C t ast , because in region- B the decaying neutrinos 
decay away when t ~ td- Furthermore, the collision terms in the Boltzmann equations for the decay radiation 
vanish for t ^ td, so that we expect that the decay radiation perturbations are well approximated by massless 
neutrinos. In this region we have the advantage that we can check this because we can calculate the CMB 
anisotropy properly. This is because region-i? the neutrinos decay out-of-equilibrium, where our Boltzmann 
hierarchy for the decay radiation, Eq. |3^, is valid. Because of this, it is possible to check the accuracy of this 
approximation. Specifically, we have checked that the calculated CMB spectrum in this region, for td -C t rec , 
is identical in the following two approaches: 1) adding a separate Boltzmann hierarchy, described by Eq. [35| for 
the decay radiation perturbations, and 2) simply increasing the effective number of massless neutrinos within a 
ACDM framework, using Eq. 

• C: m h > 3T(td), t rec <t d <t v 

In these models, the neutrinos decay out-of-equilibrium and after last scattering. The decay radiation is not 
present until after last scattering; the decays source anisotropy through the late-ISW effect. As for region-_B, 
the amount of decay radiation at decay is determined by the parameter a. But for region-C, the parameter 
degeneracy is broken, because the scale of the late-ISW feature depends on the neutrino lifetime. The CMB 
spectra for several late-decaying models are shown in Figs [| [?], |[ Note that the size of the ISW effect increases 
as rrih increases, for fixed td, and the location of the feature shifts to larger scales (smaller I) as td increases. 

We can estimate the location of the late-ISW effect by noting that it is sensitive to the scale of the sound horizon 
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FIG. 6: CMB anisotropies for late-decaying models with t d = 10 14 sec. The solid line represents the baseline ACDM model. 
The dashed line has = 10 eV, the dashed-dotted line has nij, = 31.4 eV and the dotted has = 100 eV. The decay 
radiation sources a late-ISW feature that becomes more prominent for larger masses. 
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FIG. 7: CMB anisotropies for late-decaying models with td = 10 15 sec. The solid line represents the baseline ACDM model. 
The dashed line has m/, = 10 eV, the dashed-dotted line has nih = 31.4 eV and the dotted has mj, — 100 eV. The late-ISW 
feature is shifted to larger angles relative to the td = 10 14 sec models. 



at the time the potentials decay. For neutrinos that decay out-of-equilibrium, like those in regions- £? and C, 
this time is near t — td, so that the location of the ISW induced feature is determined by the lifetime of the 
neutrino. For lifetimes shorter than the age of the universe, inhomogeneities on scales k project onto angular 
scales t ~ kr where To is the conformal time today (we assume a fiat universe). The potentials vary in time, 
and hence cause the ISW effect, most significantly at the time of decays on scales of order the sound horizon: 
k 2 h ~ 3/(4tJu;) where w = P/p is the averaged equation of state. Therefore, the bump in the spectrum is 
produced at I ~ fc s /jT — {to/tcl){^w /2>)~ 1 / 2 . If the decay occurs after matter domination but before possible 
cosmological constant domination (which occurs only at very late times), then w is determined by the decay 
radiation. Since the epoch of matter-radiation equality is near recombination for the models we are considering, 
this assumption is valid for most of region-C. Hence w ~ fi r d(*d)/3, where £l r d(td) is the fraction of critical 
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FIG. 8: CMB anisotropies for late-decaying models with t d = 10 15 sec. The solid line represents the baseline ACDM model. 
The dashed line has mj, = 10 eV, the dashed-dotted line has = 31.4 eV. The model with — 100 eV is not shown 
since this model has £l r d > 1 — Ob — = 0.22, i.e., the model is in region-i5. In these models, the late-ISW feature has 
significant power at the quadrupole, which suppresses the small angle anisotropy in this quadrupole-normalized plot. Of course, 
the normalization is allowed to vary in all subsequent analysis. 



density in decay radiation at decay. If we assume that the decay radiation never dominates the universe, then 
we can estimate £l r d(t d ) in terms of the neutrino properties: 

n r<J (t d )~1.7xlO- 8 i(^) 4/3 1 (46) 

valid for Q rd <C 1. Since we are assuming that the universe is matter dominated at decay, physical times are 
related to conformal times by r oc i 1 / 3 . If, on the other hand, £l r d(td) — T then the decay radiation dominates 
until very late times, and we have the radiation-dominated expression r oc t 1 / 2 . We can combine these results 
to obtain an approximate expression for the location of the IS W peak for region-C : 

( 1200 ha- 1 / 3 if m h < 120 /i 3 / 2 eV 
hsw^{ 1 K /i7 itm h >120h 3 / 2 eV ' (47) 

where tjj ~ 4 x 10 17 sec is the age of the universe. Entropy fluctuations, which occur when there are appreciable 
amounts of both matter and radiation, decrease the sound speed, thereby increasing lisw- The relative size of 
this effect is typically of order 20-40%. 

D: t d > tu 

In this region, the massive neutrino is effectively stable. Stable neutrinos have a long history as a dark matter 



candidate. Constraints on these models have been explored in Refs. [j49, 51 . 
E: m h < 3T(t d ), m h < 3T(t rec ) 

Here the neutrinos decay in equilibrium. Therefore, the energy density in radiation increases by SN V = 0.17 
after the neutrino becomes non-relativistic and decays away. However, since < 3T(r*), this occurs after last 
scattering, with the exact time depending on m^; the CMB anisotropy in this region are degenerate in Td- The 
small late-ISW effect that is induced is too small to be measured, even with future satellite-based experiments. 
For this reason, we will not study region- E any further. 

F: m h < T(t d ), t d > t rec 

Here, the density in either stable neutrinos or their decay radiation is enough to require £Iq > I. These models 
are extreme and suffer several problems, such as producing a universe that is too young, and so will not be 
analyzed further here. 
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VI. RESULTS 



The goal of this section is to answer two questions. 1) Is the CMB anisotropy for some decaying neutrino model 
sufficiently different from baseline ACDM so that the two models are distinguishable? 2) Given a particular decaying 
neutrino model, how well can the cosmic paramet ers, in cluding neutrino parameters, be measured? To answer question 
1) we use the distinguish abilit y framework of Sec. IV B and the Appendix, and to answer question 2) we use the Fisher 
matrix approach of Sec. IV A . 

In both cases, we adopt the following ACDM model as our baseline: £1a = 0.7, Qcdm = 0.22, Vis = 0.08, h = 0.5, 
Harrison-Zcldovich primordial spectrum (n s = 1.0), reionization optical depth t» = 0.1, and three massless species 
of neutrinos. The set of cosmic parameters allowed to vary was Xi — {^a, Qb, h, n s , t*, Q}, where Q is the overall 
normalization. To calculate the Fisher matrix onj we took two-sided derivatives for all of our cosmic parameters 
parameters as suggested in reference O], i.e., 



dC 



xi 



dX l 



Cxi{Xi + SXj) - Cxi{K - 5\j) 
28X t 



(48) 



where SXi is the numerical stepsize in the i-th cosmic parameter. All of our derivative stepsizes were taken to be 
3% of their baseline values, except for r*, whose stepsize was 0.03. We verified numerically that the derivatives were 
stable with respect to varying the stepsizes. In calculating dCxi/dflB and 8Cxi/dQ\, we allowed Qcdm to vary, to 
maintain flat universe: Qcdm = 1 — — The baseline model and its derivatives are shown in Fig. 0. 
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FIG. 9: Baseline ACDM model and its derivatives with respect to cosmic parameters. The top panel is the quadrupole- 
normalized baseline CMB spectrum: £Ib ~ 0.08, Qa = 0.7, h = 0.5, r* = 0.1, n s = 1.0. The lower panels are derivatives with 
respect to Qb, ^a, h, n a , and t„, normalized to the baseline spectrum: 1/Cxi dCxi/dXi. The derivative with respect to Q is 
not shown since dCxi/dQ oc Cxi- 

From the CMB spectrum and its derivatives, we calculated the Fisher matrix, using Eq. ^3|. To analyze a real 
experiment requires understanding details like their window functions and experimental noise. However, for future 
satellite-based experiments like MAP and Planck, the experimental uncertainty is expected to be below cosmic variance 
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Parameter 


lmax = 1000 


lmax — 2500 


Im«i = 1000 (w/ pol.) 


L«i = 2500 (w/ pol.) 




5.16 % 


2.01 % 


2.27 % 


0.700 % 




4.02 % 


1.26 % 


1.71 % 


0.391 % 


h 


3.56 % 


1.05 % 


1.46 % 


0.323 % 


n 3 


1.46 % 


0.340 % 


0.665 % 


0.194 % 


T 


0.0803 


0.0523 


0.0579 


0.0507 


Q 


5.75 % 


5.29 % 


5.37 % 


5.24 % 



TABLE I: Statistical uncertainties on cosmic parameters for the fiducial ACDM model for l max = 1000 and for l max = 2500, 
with and without polarization information. In all cases all cosmic parameters were allowed to vary simultaneously. 



for most of the angular scales they are designed to measure, and the window functions are relatively narrow. This allows 
us to characterize the experiments as cosmic variance limited to some l m ax , with the value of l ma x determined by the 
experiment. We take l ma x = 1000 for MAP and l m ax — 2500 for Planck. For both values of l m ax we consider cases with 
and without polarization information. The reason for this is that it is not certain how good polarization information 
will be. For the case that includes polarization, we assume cosmic-variance limited polarization information from a 
minimum l m i n — 200 up to the same l ma x as for the temperature data. The reason for the minimum value of I is that 
the large-scale polarization signal is small enough to be overwhelmed by the experimental noise of MAP and Planck. 
Our results are insensitive to the precise value of l m i n . The statistical uncertainties on the cosmic parameters are 
shown for MAP and Planck in Tab. VI. In this work, we take the conservative (and realistic) approach of always 
marginalizing over all cosmic parameters simultaneously. In this case, Eq. H gives the statistical uncertainty on the 
parameters. 



A. Ruling out models 

We analyzed a grid of models, consisting of 20 masses with log(m^/eV) evenly spaced from -1.0 to 1.40, and 13 
lifeti mes w ith log^d/sec) evenly spaced from 10.0 to 18.0. For each grid point, we followed the procedure given in 



Sec. rVB and the Appendix. An example of this procedure, for a late-decaying scenario, is shown in Fig. [L0[ There we 
show the ACDM and decaying neutrino spectrum, along with the best-fit perturbed ACDM model and the discrepancy 
in the fit in units of cosmic variance. From this discrepancy we calculate a confidence level for the model. The results 
for MAP and Planck are shown in Figs. |ll and |l2|. For MAP, stable neutrinos of masses greater than a couple of eV 
are distinguishable from the baseline mode , while for Planck, the sensitivity extends down to masses of several tenths 
of an eV. As the lifetime decreases and the neutrino becomes unstable, but late-decaying, the sensitivity in mass 
decreases somewhat. This is because the late-ISW signature of a late-decaying neutrino is mostly degenerate with 
rcionization. When the lifetime is short enough so that the neutrinos are decaying before last scattering, models with 
the same value of a are degenerate, and are distinguished at the same level. This is clear from a visual inspection of 
the plot. Finally, even the most optimistic case of Planck with polarization will not be able to distinguish equilibrium 
decaying models from ACDM. 

For early-decaying neutrinos, we can obtain a clearer picture by exploiting the parameter degeneracy, describing 



the models with the single variable a. Fig. 13 shows the confidence level for models as a function of a. MAP will be 
able to distinguish models with a > 10 without polarization, and a > 5 with polarization. Planck, with or without 
polarization, will distinguish any out-of-equilibrium decaying models, with a > 1. This plot confirms the result that 
models in region-A, with 8N V — 0.17 at recombination, are indistinguishable from ACDM. 



The formalism used to perform these distinguishability calculations is valid in a linear regime, where Eq. 51 holds. 
If the parameter biases become large then the formalism breaks down. Since some of the decaying neutrino models 
produce CMB anisotropy very different from the canonical ACDM, the linear approximation must break down for these 
models. However, the distinguishability contours can be believed if two facts hold. First, the linear approximation 
should hold for models that are just becoming indistinguishable, i.e., those along the contour lines in Figs, [ll] and 
|l2| . Second, models that are inside the contour must stay indistinguishable. The first point we observe to be true 
numerically The second point could break down in a couple of ways: a) the spectra start to look more like standard 
ACDM as we go inside a contour, or b) the spectra don't look like our baseline ACDM but instead look like some 
standard model with very perturbed parameters. Neither objection holds. The first is obviously false because for any 
fixed td, the decaying neutrino effects increase as we go inside the contour, increasing m^. The second objection is only 
slightly more problematic. For late decaying neutrinos, the decaying neutrino feature is a late-ISW bump at some large 
angular scale - it's pretty easy to see that this cannot be mimicked by ACDM with perturbed cosmic parameters. For 
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FIG. 10: Example of distinguishability analysis for a late-decaying neutrino scenario. Here mj, = 3.23 eV, and td = 10 sec. 
In the top panel the dashed line depicts the baseline ACDM model, with arbitrary normalization. The solid black line shows 
the decaying neutrino spectrum and the grey line shows the best-fit perturbed ACDM model. The thickness of the grey line 
represents cosmic variance. The bottom panel shows the difference between the decaying neutrino spectrum and the best-fit 
perturbed ACDM model, in units of cosmic variance. This model produces an ISW peak near I — 25, whose signature can 
clearly be seen in the bottom pane - no values of cosmic parameters in a ACDM model can reproduce such a feature. For MAP, 
without polarization, this model is ruled out at the 89.7 % level. 



early decays, the early-ISW effect is degenerate with the ratio of matter density to radiation density at last scattering. 
But in this work the only non-standard physics we are allowing is the decaying neutrino itself. This allows us to 
fix this ratio for the set of ACDM models. Other cosmic parameters affect the relative amount of radiation at last 
scattering. In particular h, and f^A, are mostly degenerate with N v jj3|. However, the degeneracy is not complete so 
that h and cannot completely mimic the early-ISW signal for these models. Since we are considering models well 
within the distinguishability contours where the early-ISW signal is large, the lack of complete degeneracy prevents 
h and Q\ from mimicking the decaying neutrino signal. 



B. Measuring neutrino parameters 

Here we are concerned with our ability to measure cosmic parameters, where the set includes quantities that specify 
the decaying neutrinos. We are primarily interested in the answers to two questions. First, what are the statistical 
uncertainties in the neutrino parameters? This goes to the goal of using the CMB as a probe of neutrino physics. 
Second, how much are the uncertainties in the non-neutrino cosmic parameters degraded by their presence? It is 
always true that adding extra parameters to the set increases or at best doesn't change the uncertainty in the existing 
parameters. If the extra parameters are orthogonal to the existing parameters in the sense that the change in the 
CMB spectrum from perturbing the new parameters cannot be mimicked by perturbing the existing parameters, 
then the degradation in the existing uncertainties is minimal. If, in the other extreme, the effect of perturbing new 
parameters can be mimicked by changing the existing parameters, the degradation is severe. Mathematically, this can 
be analyzed in terms of cross-correlations in the Fisher matrix: large cross-correlations mean degraded sensitivities. 
This degradation is one of the main arguments for pursuing the distinguishability calculations of the last section. If 
the CMB provides no evidence for decaying neutrinos, i.e., the real-universe CMB spectrum is not distinguishable 
from ACDM, then adding decaying neutrino parameters will be a hard sell. 

Since the physics behind the CMB anisotropy is different for different regions of neutrino parameter space, there is 
no one best set of neutrino parameters to add to the cosmic parameters. We will group together the early-decaying 
neutrino models, corresponding to region- A and region- £?, together, and use a as the sole neutrino parameter in this 
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FIG. 11: Decaying neutrino parameter space, showing models that are distinguishable from ACDM. The three contours represent 
distinguishability at the 90%, 99.9% and 99.9% levels for the MAP experiment. In the top panel, temperature data is considered 
alone; the bottom panel includes polarization. Models to the right of the contours are distinguishable. 



region. It is easier to compute the CMB spectra in terms of the radiation energy density N u , but a is more directly 
related to the mass and lifetime of the neutrino; the uncertainty in a can be related to the uncertainty in N v through 
the relation 

da 

5a st at = SNu,statg^r- , (49) 

where the derivative is obtained from a numerical solution to the Boltzmann equation, summarized in Fig. [}| For the 
late-decaying, but unstable neutrinos in region-C, we can just use the mass and lifetime as our neutrino parameters, 
(mh,td)- However, as the neutrinos become stable, in region-D, the inverse of the lifetime becomes a more natural 
parameter, since the baseline model corresponds to the limit v/1% — > and — > 0. Therefore, we define an 
inverse-lifetime parameter, scaled to the lifetime of the universe, 

V=%-> (50) 

td, 

and use the set (rrih, y) in region-D. Our parameter choices are summarized in Tab. 0. 

Fig. |TJ shows the relative statistical uncertainty in a, versus a, for early-decaying models. Note that 5a s tat/ot 
is a decreasing function of a. This is because the relative sensitivity to radiation energy density is roughly model- 
independent, i.e., SNv^tat/Nv is roughly constant. This means that Sa sta t/a oc a -2 / 3 , for a ^> 1. For aCl, the 
radiation density ceases to depend on a at all: Sa^t/a — > oo as a — > 0. The value of a with 8a sta t = en is interesting 
because there ACDM, with a — is ruled out at the 1-er level. For MAP, this point occurs at a ~ 25 without and 
a ~ 5 with polarization. For Planck, with or without polarization, this occurs near a = 1. Models where the data 
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FIG. 12: Same as Fig. [n], but for the Planck experiment. 



Region 


Neutrino Parameter (s) 


A 


a (degenerate, with 5N V = 0.17) 


B 


a 


C 


{m h , td) 


D 


(m h , y) 



TABLE II: Neutrino parameters to add to the set of cosmic parameters, for different regions of neutrino parameter space. 



would rule out ACDM at high significance occur for only slightly higher values of a. These values should be compared 
to the distinguishability results from the last section, where distinguishable values of a were a factor of several higher. 
This represents the advantage of including neutrino parameters in the analysis: one can rule out more models this 
way. 

This advantage comes at a price: the degradation in the ability to measure the non-neutrino parameters. First, 
consider a < 1, where the decaying neutrino models produce CMB spectra very similar to ACDM, and the uncertainties 
in the cosmic parameters reflect the degra dation that would occur if one added a (or 5N V ), and analyzed . The results 
for the limit a — > are shown in Tab VI B . Note that the relative degradation is much larger for MAP than for Planck, 
and that the degradation for certain parameters, CIb, h, n s , is quite severe. Fig. |l5|, shows the statistical uncertainties 
in the other cosmic parameters as a function of a, also normalized to the uncertainties for ACDM. The table just 
discussed is the a — > limit of the figure. The most prominent feature of the figure is general trend towards lower 
sensitivity for increasing a. 

The parameter uncertainties for several late-decaying models in region-C are shown in Tables VI B| and 
uncertainty in the neutrino parameters rrih and td increases as the lifetime increases from 10 14 to 10 1D sec. 
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FIG. 13: Level at which early-decaying models are allowed, as a function of relativity parameter a = 0.087(m h /McV) 2 (td/sec), 
for Imax = 1000 without polarization (solid line), 1000 with polarization (dashed line), 2500 without polarization (dash-dotted 
line) and 2500 with polarization (dotted line). Equilibrium-decaying neutrinos, corresponding to region- A, have a < 1, whereas 
neutrinos that decay out-of-equilibrium, region-_B, have a > 1. For two-body decays, 8N V — + 0.17 as a — > 0. An experiment 
sensitive to lmax = 2500, with or without polarization information, will be sensitive to neutrinos on the border between 
equilibrium and out-of-equilibrium decay (a ~ 1). Without polarization, an experiment sensitive only to lmax = 1000 can only 
probe very out-of-equilibrium decays (a > 100); including polarization increases the sensitivity to a ~ 5. None of the cases 
considered will be able to distinguish equilibrium-decaying models from ACDM. 




a 



FIG. 14: Using the CMB to measure a for early-decaying neutrinos, corresponding to regions-A and B. The solid lines show 
the statistical uncertainty in the parameter a as a function of a. In order of increasing sensitivity, the solid lines correspond 
to Imax = 1000 (no polarization), Imax = 1000 (with polarization), l max = 2500 (no polarization) and l max = 2500 (with 
polarization). The dot-dashed line represents the case where the 5a 3ta t = a; the dashed line shows Sa s t a t = 0.1a. For models 
below these lines, a can be measured to good relative accuracy. 



for this is that the ISW peak for the lower lifetime neutrinos occurs at higher I, where two features work to improve the 
sensitivity. First, cosmic variance is lower. Second, a given range of angular scales translates to a larger number of Vs. 
As for the other parameters, since we add two extra parameters to the analysis implies that we might expect relatively 
large uncertainties. This is observed for most parameters, especially for MAP. For some cases, the uncertainties are 
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FIG. 15: Degradation in ability to measure the non-neutrino parameters, for early-decaying models. Each panel shows the 
statistical uncertainty in a cosmic parameter as a function of a. The uncertainties are normalized to the value obtained analyzing 
ACDM without decaying neutrino parameters. The different curves in each panel correspond to MAP without polarization 
(solid), MAP with polarization (long-dash), Planck without polarization (dash-dot), and Planck with polarization (dotted). As 
a — » 0, the CMB anisotropy is close enough to ACDM so that in this limit the curves may be interpreted as the degradation 
caused by adding a as a parameter to ACDM. 



Parameter 


Imax = 1000 


Imax — 2500 


Imax = 1000 (polarization) 


Imax = 2500 (polarization) 




2.85 


1.00 


1.23 


1.21 




1.11 


1.03 


1.00 


1.65 


h 


2.21 


1.01 


1.21 


1.18 


n s 


1.52 


1.04 


1.06 


1.02 


r 


1.02 


1.00 


1.00 


1.00 


Q 


1.01 


1.00 


1.00 


1.00 
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TABLE III: Statistical uncertainties in non-neutrino parameters for a = 0, i.e., where a is added as a cosmic parameter. The 
results are shown normalized to the case where the data is analyzed without neutrino parameters. Therefore, the numbers 
represent the degradation in sensitivity from including non-neutrino parameters. Results are shown for Imax = 1000, and 2500, 
with and without polarization. 



actually less than for ACDM, which appears to violate the requirement that adding cosmic parameters decreases the 
sensitivity in the other parameters. However, here it doesn't make sense to talk about degradation, since the CMB 
spectra for these models is very different from (these models are all distin guis hable - see Figs. and H). 

The results for several almost-stable scenarios are shown in Tables VI B and VI B . The relative degradations in the 
non-neutrino parameters become worse as y increases because there the CMB anisotropy starts to be affected by the 
decay products. For very low values of y, the uncertainty degradations are the same as the case for stable neutrinos, 
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(2.53) 


3.28 
(2.24) 
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(1.02) 


80.5 
(1.39) 
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TABLE IV: Using the MAP experiment to measure mh and td for late-decaying neutrinos. The statistical uncertainties on the 
cosmic parameters, 5Xi/Xi, in percent, are shown for several models. The number in parenthesis is the ratio of the uncertainty 
to the uncertainty for ACDM. For each model the top row of data is for temperature data only and the bottom row includes 
polarization. 

with m,h as the sole additional cosmic parameter. A prominent feature of the data here is that for both MAP and 
Planck, Sy 3> y for y <C 1; it is impossible to use the CMB to probe neutrino decays in the almost-stable limit. This 
is because the late-ISW feature is imprinted at very low values of I where the cosmic variance is high and where there 
are few Vs to measure. 



VII. SUMMARY 

The goal of this work was a study of using anisotropy in the CMB to constrain the physics of neutrinos that decay 
into non-interacting daughter products. We presented the formalism required to compute the CMB anisotropy spectra 
in these models. This required calculating the energy densities and the perturbations in the decaying neutrino and 
its decay products, and incorporating this physics into the CMBFAST code [47fl . We divided the decaying neutrino 
parameter space into regions, delineated by significant physical scales, and discussed the physics behind the CMB 
spectra in each region. An enhanced early or late integrated-ISW effect is the main effect for most of the neutrino 
parameter space. 

We then developed analytic methods, valid in the linear regime, to determine when a model is distinguishable from 
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TABLE V: Same as last table, but for Planck. 



some canonical model like ACDM. With temperature data alone MAP can distinguish stable neutrino models from 
ACDM if the neutrino mass mh > 2 eV. Adding polarization data, m/, > 0.5 cV is distinguishable. Planck can 
distinguish mh > 0.5 cV with temperature alone, and mj, > 0.25 eV with polarization. MAP without polarization 
can distinguish out-of-equilibrium, early-decaying models as long as (m^/MeV) 2 td/sec > 230, and with polarization 
if (m h /MeV) 2 t d /sec > 150. For Planck without pol arization, mo dels with (m h /McV) 2 t d / sec > 9 are distinguishable, 
and with polarization if (m/j/MeV) 2 td/sec > 6. Models in which neutrinos decay in equilibrium are indistinguishable 
from ACDM. Late-decaying models (10 13 sec < < 4 x 10 17 sec) are distinguishable from ACDM if mj, > 5 eV for 
MAP and m h > 2 cV for Planck. 

Next, we studied the use of future CMB satellite data to measure cosmic parameters, including neutrino properties. 
The sensitivity to neutrino parameters depends strongly on the parameters themselves. We found that including 
neutrino parameters in a model significantly degrades the sensitivity to Qb, h, and n s , and that the degradation is 
worse for MAP than Planck. For models whose CMB spectra are not close to ACDM, the situation is less simple, 
but the sensitivities to cosmic parameters are usually less than for the canonical case. For early-decaying models, the 
sensitivities to most non-neutrino parameters decreases as a increases. In addition, we calculated the set of models 
(for early-decaying neutrinos, for now), where the statistical uncertainty in the neutrino parameters is low enough 
relative to the parameters themselves, to count as a detection of decaying neutrinos. For early-decaying neutrinos, 
MAP with can achieve this if a > 10 with temperature information alone, and if a > 3 with polarization data. The 
equivalent sensitivities for Planck are for a > 1 with temperature information alone, and a > 0.8 with polarization 
data 
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TABLE VI: Using the CMB to measure and y for nearly stable neutrinos, for MAP. The statistical uncertainties on the 
cosmic parameters, SXi/Xi, in percent, are shown for several models. The number in parenthesis is the ratio of the uncertainty 
to the uncertainty for ACDM. For each model the top row of data is for temperature data only; the bottom row includes 
polarization. 
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= 0.1 


4.64 


159. 


3.22 
(1.60) 


3.15 
(2.51) 


1.69 
(1.60) 


0.550 
(1.62) 


66.4 
(1.27) 


5.73 
(1.08) 








2.94 


62.8 


1.12 
(1.60) 


1.19 
(3.04) 


0.576 
(1.78) 


0.228 
(1.17) 


49.6 
(0.978) 


5.27 
(1.01) 




= 3.16 eV, y = 


= 1.0 


5.74 


51.0 


4.60 
(2.29) 


4.37 
(3.48) 


2.38 
(2.25) 


0.686 
(2.02) 


83.7 
(1.60) 


5.46 
(1.03) 








3.67 


20.6 


1.40 
(2.01) 


1.50 
(3.82) 


0.727 
(2.25) 


0.389 
(2.00) 


51.5 
(1.01) 


5.03 
(0.959) 



TABLE VII: Same as Tab. 



VIB| , but for Planck. 
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Although presented in the context of decaying neutrino cosmologies, the techniques developed here could easily be 
extended to more generic scenarios involving decaying particles which decay into sterile daughter products. The main 
difference in the calculation would be in determining the particle's initial abundance (the relativistic decoupling of the 
decaying neutrino simplifies the calculation in this case). Given this, the equations for the evolution of the densities 
and perturbations would be the same as for decaying neutrinos. 

In conclusion, future CMB observations promise to provide a powerful probe of neutrino physics, over a wide range 
of parameter space not easily accessed by other means. A couple of caveats are in order. First, this investigation 
was preliminary in nature. Cosmic variance limited data is a best case scenario; real-world issues like foreground 
subtraction will complicate the actual data analysis. Hopefully, the data from MAP and Planck will approach this 
ideal. Second, the real world CMB anisotropy might look nothing like any variant of CDM, with or without decaying 
neutrinos. In this case, of course, the analysis presented here would no longer be valid; one would first have to 
understand the background cosmology before going on to study the impact of decaying neutrinos. 
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DISTINGUISHABILITY OF MODELS 



This appendix deals with the question of how to determine whether or not the CMB spectra from massive decaying 
neutrino models are distinguishable from some baseline model like ACDM, which does not contain decaying neutrinos. 
Here we will restrict ourselves to the case where the decaying neutrino model produces CMB spectra that are only 
slightly different from the baseline. 

We consider the following scenario. The universe actually contains decaying neutrinos, but the experimental data is 
analyzed without considering this possibility: the set of cosmic parameters does not include vtih or td- As a result, two 
things can happen. One, the cosmic parameters measured will in general be unequal to the true cosmic parameters, 
i.e., the results will be biased. Two, the best-fit spectra may be a poor fit. If, for example, the presence of the decaying 
neutrinos changed the spectrum in exactly the same way as adding a little extra baryon density, then the measured 
baryon density would be biased, but the best fit model would fit very well. It would be impossible to disentangle the 
decaying neutrino signature from the data. We will call a model distinguishable if the best fit model is a poor one. 

To be more quantitative, we need to work through how one measures the cosmic parameters from the data. Start 
with some definitions: let {A^} be the set of cosmic parameters considered. Here, i = 1 ... N, with N the total 
number cosmic parameters. As mentioned, this set does not include nih or td- Let {A^} be the set of true cosmic 
parameters, and {AJ} be measured cosmic parameters. Finally, let {<5Ai} be the parameter biases induced by the 
decaying neutrino's, i.e., A^ = A, + SXi- The measured cosmic parameters are determined by minimizing a \ 2 statistic 
that is a function of {Aj}, given by Eqn. |3^. Here we will assume that the experimental uncertainties are just cosmic 
variance up to some maximum value of I = l max , so that the covariance matrix is that given in Eqns. || and g|. 



We know that the solution for the measured cosmic parameters with no decaying neutrinos and no noise is {A,} = 
{Xi}. Now assume that the parameter biases, SXi, are small enough so that the following holds: 

c^ orv {{K}) * c^ orv {{K}) + ^^-^ ■ ( 51 ) 

If the experiment measures the temperature anisotropy only, then X = T. With polarization information, X = T,P, C . 
This equation quantifies the statement that the anisotropy is only slightly different than for ACDM. If this holds we 
can solve for the parameter biases: 

SX k = (a,,)- 1 £ £ Syi , (52) 



l XY 



where Sxi is the "signal" : 



a s-idata ^theory fr?o\ 

bxi = C xl - C xl , (53) 
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and ctjf. is the Fisher matrix, given in Eqn. [43]. Finally, the best-fit % 2 is given by 

(nf-ttheory \ ( ^theory \ 

3x1 SXj dx-) v ™ [ Syi 5x > dk-) (54) 

To determine whether or not a model is distinguishable, we will need to understand the statistical properties of 
(Xmin)i considered as an ensemble over different realizations of cosmic variance "noise" |57|. To develop the formalism 
for estimating the contributions to Xmin from noise and signal, we break the data spectrum Cf ata into two components: 

C data = C data,0 + ^ ^ (55) 

where C^ a is the decaying neutrino spectrum without noise, and Nxi is the noise. In a perfect experiment, Nxi 
is cosmic variance. The signal Sxi breaks up similarly: Sxi — S X i + Nxi- Then the best-fit x 2 can be written 



/ aritheory \ / p.^theory \ 

^[Sxi-SXj^^ + NxtW-^h^-S^^^ + NyA , (56) 



; xy 



where SXj is the CP bias for noiseless data, i.e., without cosmic variance. 

Consider an ensemble of experiments for a given decaying neutrino model. Each experiment will have the same 
noiseless signal, but the noise, a random variable, will be different each time. Therefore, the value of Xmin wu ^ a ^ so 
vary. Associated with each value of Xmin 1S some probability that the ACDM model is allowed, denoted a. This 
probability is just the 1 minus the cumulative distribution function, £, for the % 2 distribution with l ma x — 1 degrees 
of freedom, evaluated at x 2 



\ nun ' 



otXmmJ^i-^-Umm)- (57) 



Then the confidence level for ACDM, denoted C, can be expressed as a convolution of a with the probability P(Xmin) 
of obtaining different values of Xmim 



C = 1 - 



/ d xl ml P (X 2 min) [1 -£(l m ax - 1, Xmin)] ' ( 58 ) 
JO 



To proceed further, we need to understand the shape of P(Xmin)i which is determined by the distribution of A/. 

We will treat the A/'s as Gaussian random variables with zero mean and variance determined by cosmic variance. 
In this limit P(Xmin) ^ s Gaussian too. However, the A/'s are not really Gaussianly distributed. A more realistic 
treatment [|5~H reveals that their distribution is closer to log-normal, with large high-A; tails. The disagreement is 
greater for low values of I; for high values, say with I > 50, the distribution is approximately normal. There are a 
couple of reasons why it is acceptable to approximate their distributions as normal. First, most of the statistical weight 
in distinguishing models comes from high values of I, because cosmic variance is smaller there and for experiments 
we will be considering, with l max ~ 1000, there are just more values of I that are large than small. Second, many 
different distributions, one for each A;, collectively determine the statistical properties of the Xmin distribution, and 
as the number of contributions becomes large, P(x m in) wm tend towards a Gaussian. In this case P(x m in) ^ s specified 
completely by its mean (or expectation value) \X m in)i an< ^ ^s and variance 

4 = ((Xmin) ) ~ (Xmm) 2 ■ (59) 



First, the mean. Expanding the quadratic in Eqn. |5^, we will have terms proportional to Nxi and NxiNyi and 
terms independent of Nxi- The expectation of the linear term is zero, since {Nxi} = 0, as Nxi is a Gaussian random 
variable with mean 0. For the quadratic term, (NxiNyi) = 2C xl /(2l + 1)Sxy, where Sxy is the discrete delta 
function. Therefore, we find for the expectation value of Xmim 

(n^theory \ / n^theory \ n^itheory 

S°xi ~ SX, d -^— V~ Yl S Yl - SX, 9 -^L— + E E V xx • (60) 

3 / \ 3 J i x 

This expression simplifies in certain cases. Namely, if temperature and polarization data can be considered uncor- 
rected then X ~T,P and V XY is a diagonal matrix, with Vxx ~ ftl + l)/2CjL. Then the second term on the right 
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hand side of the last equation is just equal to the number of terms in the sum, 2{l max — 1), and the expectation value 
becomes 



(x 2 ) = 2(i -ii + V V HCLLi) 

I X=T,P Xl \ 



theory \ 2 

i--^, (61) 



The variance is an unholy mess. The second term on the right hand side is the square of the mean. The first term 
looks like the following: 



k (xLJ ) = ( E E ( N *i + D ^ V xyi ( N n + ( N wm + D Wm ) V WZm (N Zm + D Zm ) ) (62) 

where 



lm WXYZ 



QQtheory 



D xl = S%-8X j U -^- (63) 



doesn't depend on Nxi- 

If the temperature and polarization are uncorrelated, this equation simplifies considerably. The sum inside the 
brackets will contain different powers of Nxi and Nx m , the objects whose expectation values are non-trivial. Note 
that if the power of either TV; or N m is odd, then that term's expectation value will vanish. In addition, terms that 
involve only Nf or have already been discussed. This allows the expression to be greatly simplified: 

:xlj 2 ) = (e e 2 ^ d A +4(u.-d(e e 1^a 2 m 

y l X=T,P XI J y i x=T,P Xl 

l,m X,Y=T,P Xl Yl 

The last term on the right hand side can be evaluated by noting the following identities: {N xl N xl ) — 3(2C xl /(2l+l)) 2 , 
and (N xl N Ym ) = (2C xl /(2l + \))(2C Ym /{2l + 1)) if X ^ Y or I ± m. Using these identities, we find a simple formula 
for the variance: 

4= 4 {Imax ~ 1) ■ (65) 

Note that the formula depends only on the number of degrees of freedom and not on (Xmin)- Finally, we can express 
Eqn. in terms of the probability distribution for Xmin> 



C =1- I d\ j= exp 

y/2ira x 



-(x 2 -(x 2 )Y 



(66) 
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